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Abstract. 

We perform Monte Carlo simulations of the hard-sphere lattice gas on the 
body-centred cubic lattice with nearest neighbour exclusion. We get the crit- 
ical exponents, (3/v — 0.311(8) and j/u = 2.38(2), where /3, 7, and v are 
the critical exponents of the staggered density, the staggered compressibil- 
ity, and the correlation length, respectively. The values of the hard-sphere 
lattice gas on the simple cubic lattice agree with them but those of the three- 
dimensional Ising model do not. This supports that the hard-sphere lattice 
gas does not fall into the Ising universality class in three dimensions. 



In this letter we study the hard-sphere lattice gas whose atoms inter- 
act with infinite repulsion of nearest neighbour pairs. The grand partition 
function is 

E v {z)=Y,z N Z v {N\ (1) 



A' 



where z is an activity and Zy(N) is the number of configurations in which 
there are N atoms in the lattice of V sites. There are many studies on: 
the square lattice § |, |, |, §, §, ft §, g [IT 
lattice i, |, I [10], |TJ, |TJ, [16 



12, 13 



17 



the triangular 
the simple cubic 



the honeycomb lattice 
lattice |, |, |, [IJ, |TJ, [Lj, the body-centred cubic lattice [|, g, [TJ, 0, and 
the face-centred cubic lattice 0. 

On the simple cubic lattice the author has estimated the critical expo- 
nents, (3/u = 0.313(9) and j/u = 2.37(2), where (3, 7, and v are the critical 
exponents of the staggered density, the staggered compressibility, and the 
correlation length, respectively ||18|| . The corresponding values of the Ising 
model are (3/u = 0.518(7) and 7/1/ = 1.9828(57) |ZT). It does not seem that 
the hard-sphere lattice gas falls into the Ising universality class in three di- 
mensions. The purpose of this letter is to obtain another piece of evidence 
of that. We estimate the critical exponents of the hard-sphere lattice gas on 
the body-centred cubic lattice. 



We carry out Monte Carlo simulations [22, ^3] of the hard-sphere lattice 
gas (|ID on the body-centred cubic lattice of V sites, where V = 2xLxLxL (L 
— 2 x n, n — 2, 3, . . . , 30), under fully periodic boundary conditions JTTJ, [TJ]. 
We measure the staggered density, 



m' 



and the staggered compressibility, 



X 



V 



V(R 2 )/A, 



(2) 



(3) 



where R = 2 (Na — Nb)/V and A^a (-Wb) is the number of the atoms in the 
A (B)-sublattice. Each run is divided into ten or twelve blocks. (• • •) is an 
expectation in a block and 7TT is one over blocks. All the simulations are done 
at the critical activity, z c = 0.7223, |20| over 12 x 10 5 Monte Carlo steps per 
site (MCS/site) or 10 x 10 5 MCS/site after discarding 5 x 10 4 MCS/site to 
attain equilibrium. We have checked that simulations from the ground state 
configuration (The atoms occupy all the sites of one sublattice and the other 



1 



is vacant.) and no atom one gave consistent results. The pseudorandom 
numbers are generated by the Tausworthe method [p4 , f25fl . 



We estimate a critical exponent and an amplitude by using the finite-size 
. For a physical quantitiy, O, we use the nonlinear chi- 
with a function of L, O(L) = AL P , where p and A are 

V 



scaling |§ [27], |8 
square fitting 



fitting parameters, p = —fijv when O 



m 1 



p = 7/1/ when O = x • We 



calculate the chi-square per degrees of freedom, x 2 /DOF, and the goodness 
|9 



Q, 



for the data set of the sizes: 



+ 2, 



2, and 



are selected so that the difference between 



of fit, 

L max . The values of L min and L max 
X 2 /DOF and 1 is the smallest. 

We get the results as follows. For m)' defined by (H), f3/v = 0.311(8), A = 
0.38(1), x 2 /DOF = 0.94, and Q = 0.50. For X V defined by (|), j/u = 2.38(2), 
A = 0.037(2), rVDOF = 0.93, and Q = 0.50. L min = 26 and L max = 48 
for rw' and x ■ Figure 1 shows the size dependence of rrv' at z = z c . The 
solid line indicates 0.38L- - 311 . Fi gure 2 shows the size dependence of x^' & t 
z = z c . The solid line indicates 0.037L 2 38 . 

The values of the critical exponents are consistent with those of the hard- 
sphere lattice gas on the simple cubic lattice, (3/v — 0.313(9) and 7/1/ = 
2.37(2) ||19|| . They are different from the corresponding values of the Ising 
model, p/u = 0.518(7) and 7/// = 1.9828(57) 

In conclusion, we have obtained strong evidence that the hard-sphere 
lattice gas does not fall into the Ising universality class in three dimensions. 
This is an interesting result since the models on the square lattice are in the 



same universality class [11, 12, 13 
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Figure captions 



Figure 1 Size dependence of the staggered density, rrtf', defined by @ at 
z = z c . The solid line indicates 0.38L -0 ' 311 . Errors are less than the 
symbol size. 

Figure 2 Size dependence of the staggered compressibility, x > defined by (|||) 
at z = z c . The solid line indicates 0.037L 2 ' 38 . Errors are less than the 
symbol size. 
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Figure 2 
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